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ABSTRACT 


A conversion factor is used to convert Geiger (GM) 
tube count rates or ion chamber currents into units of 
the incident X-ray energy flux in a' specified passband. 

A method is described which varies the passband to opti- 
mize these conversion factors such that they are rela- 
tively independent of the spectrxim of the incident pho- 
tons. This method was applied to GM tubes flown on Ex- 
plorers 33 and 35 and Mariner 5 and to ion chambers flown 
on OSO 3 and OGO 4. Revised conversion factors and pass- 
bands are presented^ and the resulting absolute solar X-ray 
fluxes based on these are shown to improve the agreement 
between the various experiments. Calculations have shown 
that, although the GM tubes on Explorer 33 viewed the sun 
off-axis, the, effective passband did not change appreciably, 
and the simple normalization of the count rates to the count 
rates of a similar GM tube on Explorer 35 was justified. It 
is also shown that for the type 6213 GM tube, the passband 
is effectively 3 to 12 A although gas pressure and window 
thickness variations . may give rise to numerical conversion 
factor differences of a factor of 6. The programs used in 
these calculations are given in the appendixes. 





SECTION 1. APPROACH TO THE PROBLEM 


Geiger (GM) tubes and ion chambers have been used to measure the 
broadband X-ray flux from the sun. However, in using these devices it 
is necessary to convert the instrumental output from GM tiibe count rates 
or ion chamber currents into scientific .units in the form of incident 
energy fltixes. Usually, the instrumental output is multiplied by a con- 
version factor to effect this change. Because b£ different approaches 
to solving the problem of detemining the conversion factor, data ap- 
pearing in the literature are not always consistent. 

One approach to finding the conversion factor has been to assume 
both a spectral shape for the incident spectnmi and a passband for the 
detector and then to calculate the conversion factor using these assump- 
tions (_e.g. , Kreplin, 1961). But the conversion factor/passband combi- 
nation so chosen may be very sensitive to the shape of the input spec- 
trum, which in turn is hot measured by the detector and is known to 
vary. 


In another approach, developed by Van Allen (1967), a number of 
input spectra are chosen and the passband limits are varied until lim- 
its are found which minimize the spectral dependence of the conversion 
factor. This approach, which minimizes the dependence of the resulting 
data upon a quantity which cannot be measured by the instrument, was 
chosen for further development. 

The efficiency of a GM tube or ion chamber, e, equals the .proba- 
bility that a photon of wavelength X will be transmitted through the 
window -and .absorbed by the filler gas . The efficiency is determined' 
theoretically by the product of the window transmission, T^(X) , and 
the abso 2 T)tion of the filler gas, Ag(X). For a nonreflecting gas, the 
sum of the gas transmission, Tg(X), and the absorption, Ag(X), is unity. 
Thus, 

e(X) = T/X) [l - TgCX)] (1) 

where 

T(X) = exp - (ij/p)px, 

Cu/p) = the mass absorption coefficient, and 
px = the thickness of the window, w, or 
the gas, g, in g cm“^. 

The mass absorption coefficients used in this study were taken from Henke 
and Elgin (1970) and McMaster et al. (1969) aitd^were interpolated where 
necessary assuming a power law relationship between the mass absorption 
coefficient and the wavelength. (Appendix A briefly describes the pro- 
gramming language used, APL on the IBM system 360; Appendix B describes 
the program.) 



The conversion factor, y, is defined as the' ratio of the energy 
flux within a specified passband to either the GM tube' count rate or 
the ion chamber current. The energy flux is the integral, between 
specified limits (e.g., a and b) , of the differential- energy flux, 
dE/dX. The GM tube count rate equals the window area. A, times the 
integral from zero to infinity of the differential photon flux, dN/dX, 
times the probability that a photon is absorbed, e. The conversion 
factor, y, resulting for a GM tube is given by 
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CdE/dX) dX 


Tr. * 
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CdN/dX) ♦ s ; dX 


( 2 ) 


The ion chamber current equals the window area. A, times the gain of 
the gas, G, which is the number of ion-electron pairs' formed per erg 
of absorbed energy, times the charge of an electron, e, times the en- 
ergy per absorbed photon, hc/X, times the integral from zero to infin- 
ity of the differential photon flux dN/dX, times the probability of 
absorption e. The conversion factor, y, resulting for an ion chamber 
is ^ven by 


Yt 
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Chc/X) (dN/dX)-s,.dX 


C3) 


Physically, the differential energy flux is equal to the product of the 
differential photon flux and the energy of each photon, or.’ 

dB/dX = Chc/XD CdN/dX) ■ . C4) 

Applying equation (4D , equations (2) and ^3) can be cast entirely 
in terms of either the differential photon flux dN/dX or the differ- 
ential energy flux dE/dX. The conversion factor y is a function of the' 
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shape of the spectrum used and not of the magnitude of the spectrum 
since constant factors cancel out. The difference in spectral re- 
sponse between using a tube as a GM tube or as an ion chamber appears 
as a factor of 1/X in the integral in the denominator since equation 
(3) is equivalent to equation ( 2 ) if the .efficiency In equation (2) is' 
replaced by the product of the efficiency, the photon energy hc/X, the 
gas gain G, and the electronic charge e. 

A computer program has been written which uses the Van Allen cri- 
terion (i.e., spectral independence of the conversion factor) to deter- 
mine the best conversion factor/passband combination for a given detec- 
tor. This program is described in Appendix C. The input spectra chosen 
are free-free and blackbody spectra at nominal temperatures of 1-, 2-, 

5-, 10-, 20-, 50-, and 100-million °K. When con^ared to typical solar 
spectra, these spectra vary from those which are too soft through those 
which are too hard. Two flat spectra, dE/dv = 1 and dE/dX = 1, are used 
as controls. Given the efficiency of the tiibe, the window area of the 
tube, the gas gain if applicable, and the estimated limits, the program 
calculates 16 conversion factors, y's (one for each spectrum), their 
mean, the standard deviation, a, and the relative error, o/mean. The 
program varies the limits and finds the passband which gives the mini- 
mum relative error. Integrals are calculated by summing in 0.25-A steps; 
the denominator is integrated from 0,25 to 20 A, 

While the spectra used are not ideal models of the solar spectrum, 
it is expected that the conversion factor/passband combinations deter- 
mined through their use will not vary significantly from those obtained 
using actual solar spectra.' The nominal uncertainty in the passbands is 
+0.25 A, and the relative error is given for each conversion factor. The 
effect of using these conversion factor/passband combinations on actual 
data from different experiments is shown in the next section. 

Two further questions need to be asked, particularly in the case of 
the University of Iowa experiments flown on Explorers 33 and 35. These 
experiments xosed the same type of GM tube (EON 6213) , but none was ac- 
curately calibrated for soft X-ray use prior to launch. All of the data 
were normalized to the data from one GM tube on Explorer 33. Although 
these GM tubes viewed the sun off-axis, the data were again normalized 
using empirically determined angular co3rrection factors. The validity 
of these procedures is dependent upon the passband of the GM tubes re- 
maining relatively constant regardless of the known window thickness and 
gas pressure variations and regardless of the angle between the sun and 
the axis of the GM tube. 
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The program explained in Appendix D has been developed to deter- 
mine the effects of window thickness and gas pressure variations. In 
this program, a passband is assumed and conversion factors and their 
corresponding relative errors are calculated. The results of xising 
this program are discussed in Section 2. 

The program e:q)lained in Appendix E has been developed to deter- 
mine the effect of viewing a source off-axis. In this program, the 
geometry of the aperture covering the GM tube is used for input, and 
the relative efficiency and area are calcula.ted as a function of the 
off-axis angle. From these variables, the relative off-axis response 
is calculated. Again, the program validates the procedures used in 
the reduction of Explorer 33 and Explorer 35 data. 
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SECTION 2 . APPLICATIONS 


Data from six detectors on four spacecraft were used in this study. 
The detectors were the GM tubes' flown on Explorer 35 and Mariner 5 which 
were used to observe solar X rays (University' of Iowa_ experiments ; J. A. 
Van Allen, Principal Investigator), the ion chamber e: 2 q)eriment flown on 
OSO 3 (University of Michigan experiment; Prof. R, Teske, Principal In- 
vestigator) , and three ion chambers flown on OGO 4 p^aval Research Lab- 
oratory experiment; R. Kreplin, Principal Investigator) . Although data 
from Explorer 33 were not used in 'this study, the GM tubes it carried 
were the same type as those flown on Explorer 35. All of these detec- 
tors measured the broadband solar X-ray flux .in different passbands. 

First, the optimal pas sb and/ conversion factor combinations were, 
calculated. These are listed in Table 1. Applying the revised pass- 
band/conversion factor combinations yields the following changes . 

• The Explorer 35 fluxes should be multiplied. by a factor, of 1.37 
and the passband decreased to 2.75 to 12 A. 

• The Mariner 5 fluxes should be multiplied by a factor of 1.85 
and the passband increased to 2 to 9.75 A. 

• The OSO 3 flux values should be multiplied by a factor of 2.5 
and the passband increased to 1.5 to 13 A. 

• The OGO 4 8- to 20 -A fluxes should be multiplied by a factor of 
1,05 and the passband shifted to 2 to 14 A. 

• The OGO 4 1- to 8 -A fluxes should be multiplied by a factor of 
0.28 and the passband decreased to 0.75 to 7 A. 

• The OGO 4 0.5- to 3-A fluxes should be multiplied by a factor of 
0.83 and the passband increased to 0.25 to 3 A. 

The reason for the drastic change in the passbands of the OSO 3 8- to 
12-A detector and the OGO 4 8- to 20 -A detector is that the peak in the 
efficiency at wavelengths shorter than the K- edge of the- aluminum in 
the window (7.94 A) had been neglected previously due to the soft spec- 
trum (2 million °K gray body) used. This spectrum contained relatively 
few photons at the shorter wavelengths . . 

Both new and old conversion factor/passband combinations were 
applied to data obtained during times of low and high solar activity 
and during a flare. (See Table 2.) Using the old conversion factor/ 
passband combinations, numerous inconsistencies occurred. Typically, 
the flux observation by OSO 3 in the 8- to 12-A passband was about 
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Table 1, Old and Revised Conversion Factor/Passband Combinations 


Experimenter 

Spacecraft 

Old 

Passband (A) 

Old 

Conversion Factor 

Revised 
Passband (A) 

Revised 

Conversion Factor 

• Revision 
Factor, 

Univ. of Iowa 
Explorers 33 ,35 

2-12 

1.8x1 0'^^ ergs cm"^ cf’ 
3 . 5x1 0 ^ ergs cm '2 ct‘^ 

2.75-12 

(2.47±.29)xl0'^ 

1.37 

Mariner 5 

2-9 

2-9.75 

(6.46i.51)xl0 

1.85 

Univ. of Mich. 
OSO 3 

8-12 

7,71x10® ergs cm'^amp'^^^ 

1.5-13 

(l,93±,43)xl0’ 

2.5 

Naval Research 
Laboratory (NRL) 
OGO 4 

8-20 

5.25x10™ ergs cm'^amp'^ 

2-14 

(5.5±1.25)xl0™ 

1.05 

OGO 4 

1-8 

5.7x10® ergs cm'^amp’’ 
1.48x10® ergs cm’^amp’* 

0.75-7 

(1.59±.6)x10® 

0.28 (5) 

OGO 4 

0.5-3- 

0.25-3 

(1.23 £.17)xl0® 

0.83 


(1) (Old flux value) =(revised flux value) 

(2) Uncertainty of ^50% 

(3) Uncertainty of -30% 

(4) Spectral temperatures 5, 7.5, 10, 25, 50, 75, & 100 x 10^ “K; 
all others, 1, 2, 5, 10, 20, 50, & 100 x 10*' *K 

(5) Using the temperature in (4), the relative error decreases from .38 to 
.057, the passband decreases to 1- to 5.75- A, the conversion factor 
becomes (1 . 145 i ,065)xl0®, the revision factor becomes 0.20. 




Table 2. Typical Observed X-Ray Fluxes 









twice the flxix observed by Explorer 35 in the 2- to 12-A passband, im- 
plying the nonphysical condition of a negative solar X-ray flttx between 
2 and 8 A. Further, the OGO 4 1- to 8-A flux nearly equalled the' Mari- 
ner 5 2- to 9 -A flux during hi^ activity (i.e., the flux came to within 
the. uncertainty in the Mariner 5 flux) . During a flare, the OGO 4 1- 
to 8-A flux exceeded both the 2- to 9-A flux obseirved by Mariner 5 and 
the 2- to 12-A flux observed by Ejq>lorer 35. Employing the revised con- 
version factor/passband combinations eliminates these inconsistencies. 
These integral flxixes are illustrated in Figures 1, 2, and 3. 

Since the solar spectrum in this region decreases strongly with de- 
creasing wavelength, these flux measurements can be combined to yield 
crude integral spectra. It must be assumed that the major contribution 
of the incident flux to the detector response occurs at the long wave- 
length edge of the detector passband. Such integral spectra are shown 
in Figure 4. The error bars in Figure 4 indicate the relative error 
calculated. Using the ratios of GM tube count rates or ion chamber cur- 
rents, crude spectral temperatures can be determined. (See Appendix F.) 
Temperatures obtained using this method range from 4 to 8 million °K for 
wavelengths less than 9 A to about 1 million ®K for longer wavelenths. 

Having shown that the data bases can be made consistent, the integ- 
rity of the Explorer 33 and Ejfplorer 35 data bases is questioned next. 
Figure 5 shows the conversion factor/passband combinations derived for 
a number of different window thicknesses and gas pressures. The pass- 
bands seem to cluster about a nominal passband of 3 to 12 A, Setting 
a nominal passband of 3 to 12 A increases the worst case of relative 
error to about 0.4', roughly the value estimated by Van Allen. Given 
this relative error, the normalization of all the data to the data from 
one GM tube on Explorer 33 seems justified. 

The question concerning the off-axis response was answered using 
the program given in Appendix E. The angular correction factor for GM2 
and GM3 on Explorer 33 was matched against the inflight value every 2.5" 
from 0“ to 25" . The geometry was varied slightly until the match was 
within 10 percent at every pmnt. The effective efficiency determined 
every 2.5° was used as input to the program, which determined the best 
conversion factor/passband combination, and the optimal passband was 
found. The results are shown in Figure 6. Note that the passband 
shifted only 0.25 A at the short wavelength edge and that a shift of 
this level is at the level of integration noise since the integrations 
were summations in 0.25 -A' steps. An alternative means of showing the 
spectral independence of the off-axis response is by taking the ratio 
of the effective efficiency obtained when viewing the tube at different 
angles to the effective efficiency obtained when viewing the tube head- 
on.. In the region of 3 to 12 A, these ratios (Figure 7) are relatively 
constant. Thus, the empirical normalization procedure is valid. 
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Figure 1. 


PASSE 

Integral Flux on Sepfembi 
(Low Solar Activity). 





INTEGRAL FLUX (ergs cm~?sec'J) 



Figure 2 . Integral Flux on August 20 , 1967, 1200 to 1300 UT 
(Moderate to High Activity, No Flares) 
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MNitVartAL TLUA um 



PASSBAND (A) 


Figure 3. Integral Flux on August 19, 1967, 0730 to 0800 UT 
(Post-Peak Decreasing Phase, 

3-12 A Flare 0557 - 0615 - 0833; 

2B Optical Flare 0530 - 0605 - 0630; 

2N Optical Flare 0730 - 0735 - 0800) 



INTEGRAL FLUX (ergs cm-^sec-^) 



Figure 4. Integral Spectra 



A -28 cm OF Hg OF NEON 
B-56 cm OF Hg OF NEON 
C-112 cm OF Hg OF NEON 


WORST CASE - 14% UNCERTAINTY 


-NOMINAL- 


-1.75A 

T75B 

-L— 1.25 A 
•1.75C 
i-1 — 1.25B 

-Li.25C 


PASSBAND (A) 


16 18 


Figure 5. Comparison of Responses of E'OIN 62J.3 and Similar GM Tubes. 
EON 6213 and- Anton 213 - 1.25 mg cm'^ Mica Windows 
LND 704 - 1.5 to 2.0 mg cm"^ Mica Windows 




CONVERSION FACTOR (ERGS/cm^-count) 



2 4 6 8 10 12 14 

PASSBAND (ANGSTROMS) 


Figure 6. The Passband and Convereion Factor as a Function of 
0.ff-Axis Angle a (Explorer 33, GM2 and GM3) 
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AREA X EFFICIENCY RELATIVE TO 0* OFF-AXIS 
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APPENDIX A. A BRIEF DESCRIPTION OF THE APL LANGUAGE • 


This appendix is intended to provide the definitions necessary to 
change the programs listed in the following appendixes into other com- 
puter languages. More complete descriptions of the APL language and 
its use are found in Palcin (1968) and Gilman and Rose (1970).. 

General Characteristics : 

Accuracy - Double precision is used except for integers, which do not 
require it. 

Order of Execution - Operations are processed from right to left, with 
operations enclosed within parentheses being performed first. 

Thus, A - B - C would be equivalent to A (B - C) . 

Classes of Operators - There are two classes of APL operators, monadic 
(right argument only) and dyadic (both right and left arguments) . 
The same symbol may be used for both a monadic operator and a 
dyadic operator, but ther-meanings will be different. For example, 
*X means log^X, while 10*X means log^^X. 

Classes of Variables - Variables may be scalars, vectors, or matrices. 
They are described as follows: 

Scalars - rank 0, no dimensions, no components. 

For example, X = 9 

Vectors - one dimension, one or more components . 

For example, X =1 36 4 5 

Matrices or Tensors -- more than one dimension, one or more 
components in each dimension. 

For example, X = 1 36 

2 7 4 

9 5 7 

Most operators will function with or on any variable. If 
two or more variables are involved, their ranks must match 
(both dimensions and components) . The exception to this 
requirement is the scalar, which may be used with vectors 
or matrices Oe.g., in multiplication, addition, etc.). In 
this case, the scalar is automatically matched with the vec- 
tor or matrix. 


A-1 



Table A-1 describes the operators used in the following programs. 
The variables used are A and B, which may be scalars, vectors (with 
components Aj and Bj) , or matrices Cwith components Ajj and Bjj, 
although matrices may have more than two dimensions) . Components 
are identified as follows; 

Scalar - no components 
Vector - A[I] 

Matrix - A[I;J] 

Tensor - A[I;J;K;1] 



Table A-l . APL Operator Definitions {Sheet 1 of 3) 


OPERATOR 

C^~B 

C^A-B 

C^A>iB 

C->rAiB 

C*-*B 

C-<^A*B 

C*-&B 

C*-A9B 

C^\B 

C^A\B 

C^IB 

C^AilB 

<7-6- 1 B 
C^A IS 

C-<rOB 

(7-6-IoB 

C-6-20B 

( 7 - 6-305 

<7-«-'ioS 

(7-6--20B 

(7-6-“ 30S 

C^A<B 

C^A<B 

C^A=B 

C-<rA*B 

C<-A>B 

C^A>B 


RESULT STORED AS C 

THE VALUE OF B 
TEE SUM OF A AND B 

THE VALUE OF B I-/ITH THE OPPOSITE SIGN 
THE' DIFFERENCE OF A MINUS B 

THE SIGN OF S - IF B<0, 0 IF S=0, 1 IF B>0 
A MULTIPLIED BY B 

THE RECIPROCAL OF B 
A DIVIDED BY B 

EXPONENTIAL TO TEE BASE E OF B 
A TO THE POUER B 

THE LOG. TO THE BASE E OF B 
THE LOG. TO THE BASE A OF B 

THE NEXT INTEGER >B 
THE LARGER OF A AND B 

THE NEXT INTEGER 
THE SMALLER OF A AND B 

THE ABSOLUTE VALUE OF B 

THE REMAINDER OF B DIVIDED BY A 

FI TIMES B 
SIN B. B IN RADIANS 
COS B, B IN RADIANS 
TAN S, B IN RADIANS 
SIN'l B, C IN RADIANS 
COS“’l B» C IN RADIANS 
TAN~1 By cm RADIANS 

■C-1 IF A IS LESS THAN B, IF NOT, 0=0 
0=1 IF A IS LESS THAN OR EQUAL TO B , IF NOT, 0=0 
0=1 IF A EQUALS B, IF NOT, 0=0 
0=1 IF A DOES NOT EQUAL B, IF NOT, 0=0 
0=1 IF A IS greater THAN OR EQUAL TO B, 

IF NOT, 0=0 

0=1 IF A IS GREATER THAN B, IF NOT, .0=0 
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Table A-1. APL Operator Definitions (Sheet 2 of 3) 


OPERATOR 

C’cAaB 

C * — B 
C-eiB 


C+iB 

C^A\B 

C-^^A 

C-^A^B 

C-^-pB 

C*-ApB 

C^^B 

<7+6?B 

C-*-A ® > 
C^-i-/B 

C+x/B 

t7-t*x/S 


dPSBiiyjoiv 

LOGICAL A AETt B U*B = 

LOGICAL A OR B OR BOTH U»B = 

LOGICAL A OR B BUT NOT BOTH U.B s 1»B X7/B 0»S) 
(NOT) B, INTERCHANGES I’B 47I7B O'B 

C-t- SUBSCRIPTS OF THE VECTOR B FOR AN ASCENDING 
SORT 

C-c SUBSCRIPTS OF THE VECTOR B FOR A DESCENDING 
SORT 

C* A VECTOR OF INTEGERS FROM 1 TO B 

C*- THE INDEX OF THE FIRST ELEMENT IN THE VECTOR A 

EQUAL TO THE VALUE OF B 

C<- A VECTOR OF THE ELEMENTS WHICH COMPOSED 
TEE SCALAR, MATRIX OR TENSOR A 

C*- A VECTOR COMPOSED OF THE VECTOR B APPENDED TO 
TEE END OF THE VECTOR A 

C<- THE DIMENSION OF B , A BLANK IF B IS A SCALAR, 

A SINGLE NUMBER IF B IS A VECTOR, 

A VECTOR OF NUMBERS IP B IS A MATRIX OR TENSOR 
C^ A TENSOR, MATRIX, OR VECTOR OF DIMENSION A FROM 
THE VALUES OF THE SCALAR OR VECTOR B 

ClIl-^BlN-Il, WHERE N IS THE NUMBER OF COMPONENTS 
IN B {REVERSES THE ORDER OF ELEMENTS IN B) , OR 
Clhll^BlIiN-Jl, WHERE N IS THE NUMBER OF 
COMPONENTS IN THE ROWS OF MATRIX B {REVERSES THE 
ELEMENTS IN EACH ROW OF A MATRIX) 

C’CJiJl-f-BCc7;J]. TRANSPOSES THE MATRIX ROWS/ COLUMNS 

23 ALJl^BLJin OR CLliKl^J^ All iUxBU iKl 
J J 

ClhJ:i*-All2iBl/n, GENERATES A MATRIX WITH VALUES 
FOUND IN VECTORS, OTHER OPERATORS MAI BE USED 

C^ ilCJ] OR CCX3-*-2^ 

J 

c<r -4Cr: OR ciii^J2 All wi 

' I 

Bin OR n buwi 

. J 

c^ Bcn OR cm*- n Biim 

I 
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Table A-1 . APL Operator Definitions (Sheet 3 of 3) 


OVER A TOR OVER A TIO N 

C-i'l/B C IS THE LARGEST ELEMENT IN THE VECTOR B , ' OR A 

VECTOR OF THE LARGEST ELEMENTS IN THE RONS OF 
THE MATRIX B 

C^UB C IS THE LARGEST ELEMENT IN THE VECTOR B , OR A 

VECTOR OF THE LARGEST ELEMENTS IN THE COLUMNS OF 
TEE MATRIX B 

C^l/B C IS THE SMALLEST ELEMENT IN THE 

VECTOR OR A VECTOR OF THE SMALLEST 
ELEMENTS IN THE RONS OF THE MATRIX B 

C*-\.fB C IS TEE SMALLEST ELEMENT IN THE 

vector's, or A VECTOR OF THE SMALLEST 
ELEMENTS IN THE COLUMNS OF MATRIX B 

CALLS- FOR INVUT FROM THE TERMINAL OF A 
VALUE FOR X 

U^^X 'DISPLAYS X ON THE TERMINAL 

NAME: PRECEEDS AN INSTRUCTION AND LABELS 

A LOCATION OR LINE NUMBER WITHIN A 
PROGRAM (THE NAME USED IS ARBITRARY) 

■^NAME BRANCHES THE PROGRAM TO THE 

LOCATION CALLED NAME 

■yiLOGlO/NAME IF THE LOGICAL STATEMENT IS TRUE, TEE 

PROGRAM BRANCHES TO THE LOCATION 
CALLED NAME, IF NOT TRUE THE PROGRAM ■ 
PROCEEDS TO THE NEXT LINE 

■^{LOGIC J^J-^D/NAME INCREMENTS J, TESTS THE VALUE AND 

BRANCHES TO LOCATION CALLED NAME IF 
TEE LOGIC IS TRUE, IF THE LOGICAL 
STATEMENT IS UNTRUE, THE PROGRAM 
PROCEEDS TO THE NEXT LINE 

TIME OF DAY IN 1/60 SEC, 

CENTRAL PROCESSOR UNIT (CPU) TIME 
USED THIS SESSION (1/60 SEC,) 

NUMBER OF TERMINAL USERS ON SYSTEM 


120 

I21 

123 



APPENDIX B. MASS ABSORPTION COEFFICIENTS 


The program called INTERP interpolates or extrapolates the needed 
mass absorption coefficients that are not already listed in the tables 
of Henke and Elgin (1970) or McMaster et'al. (1969). The program is 
based on the assumption that the mass absorption coefficients are re- 
lated to the wavelength through a power law. The geometry employed is 
illustrated in Figure B-1. In Regions I and III, the values are extra- 
polated while, in Region II, the values are interpolated. The program 
is flowcharted in Figure B-2 and listed at the end of this appendix. 

The operation of the program is as follows. The data are entered 
and converted to natural logarithms (i.e., the data are linearized). 

The slope of the graph of known mass absorption coefficients versus 
wavelengths is determined between each successive pair of points. The 
slopes are corrected for the effect of absorption edges, and a multi- 
plier is determined to correct for edges found in Regions I and III and 
for the extrapolation in Region IV. Using the slopes, the multipliers, 
and the known points, the required coefficients are calculated. 

In detail, the program operates as follows. 

Steps 1-4 require the entry of the known coefficients and their 
corresponding wavelengths and convert them to natural 
logarithms. The data must be in the form of vectors. 

Steps 5 § 6 order the above data 'by ascending wavelength. 

Steps 7-12 require the entry of the known edges, the corresponding 
jump ratios, and the wavelengths of the wanted coeffi- 
cients and convert them to vectors of natural logarithms . 

Steps 13 - 15 determine the number of successive pairs of points, the 
slopes (a vector) between successive pairs of points, 
and set the multiplier as a vector of zeros, one com- 
ponent corresponding to each coefficient wanted. 

Steps 16 § 17 set up a loop which cycles through each absorption edge 
and makes the necessary corrections to the slope or to 
the multiplier. Step 17 branches to the end of the loop 
if the edge lies outside the region of interest. 
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Step 18 


corrects the multiplier i£ the edge falls in either Region 
I or* Region III (dashed lines in Figure B-l) by adding 
(Region I) or subtracting (Region III) the logarithm of 
the jump ratio to coefficients at wavelengths shorter 
(Region I) or longer (Region III) than the edge. 

Step 19 branches to the end of the loop if the edge falls within 

Region I or Region III. 

Steps .20 - 22 modify the slope and the multiplier for edges falling in 
Region III and wanted points (multiplier) falling in 
Region IV. 

Step 23 increments the index J, tests it, and loops back to 

step 17 if necessary. 

Steps 24 — 26 calculate the answer using the slopes, the known coeffi- 
cients, and the multipliers and convert the answer out 
of its logarithmic form. 

Steps 27-31 convert the initial data out of its logarithmic form. 

Step 32 announces that the program is finished. 
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STEPS 1-12 


STEPS 13 - 15 


STEPS 16-17 


STEP 18 


STEP 19 


STEPS 20-22 


STEP 23 


STEPS 24-31 


STEP 32 


Figure B-2^ Flowchart for I NTERP Program 
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V INTERP 

[11 CENTER KNOWN COEFFICIENTS'.' 

[21 CKN^^CKN-frU 

[31 » ENTER CORRESPONDING WA 7ELENGTUS : ' 

[4l WKN*-eWKN^U 

[5l CKN^CKNikWKNl 

[ 6 ] WKN^WKNlkWKNl 

[71 'ENTER WAVELENGTRS OF EDGES'.' 

[81 EDG'f,®EDG'frU 

[91 'ENTER CORRESPONDING JUMP RATIOS:' 

[lOl JR-^,9JR*-D 

[111 'ENTER WAVELENGTHS OF WANTED COEFFICIENTS:' 

[12l W'*-t®W-^U 

[131 c7-f*"l + pC7Xi7 

[141 SL'(-(,CKNtl-i-iJl-CKNL\Jl)i(WKNll + xJl-WKNL i«7^1 ) 

[151 MULT-<-{pW)pO 
[161 M 

[171 SLMOD:HiEDGlJJ<WLll)+(.EDGlJl>WlpWl))/ENDMOD 
[181 MULT<rMULT+JRU:\xUCEDGZJl<WKNLlJ)x(EDGLJl>W))-iEDGtJ'} 
>WKNlpWKNl)x(.EDGlJ2<W)) 

[191 M(EDGLJl<WKNLll)i-(.EDGlJl>WKNipWKNl))/ENDMOD 
[201 I<-+/WKN<EDGLJ1 

[211 SLLlJ-^SLZn + URLJliiWKNlI+ll-WKNLll)) 

[221 MULT-^MULT-iiW>EDGlJl)^(,W<WKNlI+i:\))^JRLJl 
[231 ENDMOD :^UpEDG) >JWi-l ) /SLMOD 
[2 41 lHl=0)-i-I<-f/Wo .■S:WKN 
[251 SL'f-SL^SLLpSLl 

[261 ANS-^*(.MULT+iCKNLIJ+SLLllxiW-WKNil2) )) 

[271 W^*W 
[2 81 CKN-<r*CKN 
[291 WKN-c*WKN 
[30l JR<-*JR 
[311 EDG^*EDG 

[321 'ANSWER IS STORED AS "ANS".' 

7 
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APPENDIX C. PASSBAND DETERMINATION PROGRAM 


The program called PASSBAND determines the optimum passhand for a 
GM tube or an ion chamber following the procedure outlined in Section 1, 
•'Approach to the Problem." The program is flowcharted in Figure C-1 and 
listed at the end of this appendix. 

The . operation of the program is as follows. The GM tube efficiency 
or the ion chamber efficiency (multiplied by the gas gain, electronic 
charge, and photon energy) is entered followed by the area and the esti- 
mated limits. The temperature range- is chosen, and lines are added to 
the spectra, if appropriate. For each pair of limits, set up as a ma- 
trix with -the lower limit varying from the estimated lower limit -0.5 A 
to +0.5 A (in 0.25-A steps) versus the estimated upper limit -0.5 A to 
+0.5 A (in 0.25-A steps), the program determines the 16 conversion fac- 
tors, Y's, the mean conversion factor, the standard deviation, a, and 
the relative error, o/raean. A matrix of the relative error is derived 
(in effect, if not in fact), the minimum relative error is determined, 
and the corresponding limits are found. If the limits occur on any 
edge of the matrix, the estimated limits are jumped so that the oppo- 
site edge of the next matrix will correspond to the edge, the estimated 
limits are set equal to the limits corresponding to the minimum, the 
matrix is recalculated and displayed, and the final limits are printed 
out followed by the linearly averaged conversion factor and its corre- 
sponding standard deviation, the logarithmically averaged conversion 
factor and its standard deviation multiplier, and the 16 conversion 
factors (one for each spectrum) calculated for the optimum passband, , 

If , the minimum occurs twice in succession and has the same value, the 
program sets the estimated limits to the limits corresponding to the 
minimum and cycles through to the final results described above. 

This procedure prevents any infinite oscillation from occurring 
should the final result lie on the common edge of two adjoining matrices, 
while permitting the program to jiwip in the direction of the absolute 
minimum rather than increment one step at a time. In the rare case 
where one limit falls within 0.25 A of either 0 or 20.25 A, the final 
minimiM cannot occur in the center element of the relative error matrix. 
(This condition would require undefined values of spectra dN/dA and ef- 
ficiency.) This condition is recognized because the minimums are iden- 
tical three successive times, and the program prints out the off-centered 
matrix with the final results . 
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In detail, the program operates as follows.- 

Steps 1-4- initialize variables. 

Steps 5-8 require the entry of the GM tube efficiency (or the ion 
chamber equivalent) , multiply it by the inverse of the 
photon energy X/hc, and convert it to an integer accu- 
rate to one half part in- 10® with a corresponding nor- 
malization factor XY. 

Steps 9-12 require the entry of estimated limits (a two-element vec- 
tor) and determine the corresponding indices in a stan- 
dard vector of wavelengths . 

Steps 13 - 19 require the entry of the desired temperature range (en- 
tering the letter Y enters the values set in step 13) and 
calculate relative free-free and blackbpdy spectra for 
the temperatiires . Step 16 determines a matrix of values 
of hc/WkT, with the wavelength W in rows versus the tem- 
perature T in columns . Step 17 sets a maximum value of 
170 on independent elements in the matrix, and steps 18 
and 19 determine relative spectra (free-free spectra are 
determined in step 18, blackbody in step 19) in integer 
form accurate to one half part in 10®. 

Steps 20 - 23 require the entry of the wavelengths of lines (each, in 
turn, determining the corresponding index of the wave- 
length vector IV) and the corresponding line/ continuum' 
ratios , 

Steps 24 - 27 form a loop (step 24 sets the index as J = 1) which multi- 
plies each element in a row corresponding .to a line wave- 
length by a factor of the line/ continuum ratio + 1. Step 
27 increments 'the index J and tests its value against the 
total number of lines. 

Step 28 generates a 16 -element vector of count rates, one for each 
spectrum. The integer form of the efficiency is corrected 
by the factor (XY • lO”®). The terms like (GM +. x FF) are 
equivalent to 

I 

The last term adds •die spectra dN/dX = 1 and dN/dX = X” . 

Step 29 clears the value of GM. 
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Steps 30 - 32 check that the upper and lower limits estimated are not 
too close to each other or to either 0 or 20.25 A. Step 
30 determines the difference between the two limits, X. 

Step 31 accomplishes several functions. First, if the 
limits are too close together (X < 4) and the estimated 
lower limit 10 > 8, then the difference between 4 and X 
is subtrated from the value of 10. Second, if 10. > 2, it 
retains its current value. Finally, if 10 < 3, it is re- 
placed by a value of 3. Step 32 provides that if 10 5. 8, 
the difference between 4 and X is added to JO. If JO < 79, 
it retains its current value. If JO > 78, it is set equal 
to 78, 

35 set I and J. Step 35 clears variables. 

40 determine a 16 -element vector (GV) of energy fluxes di- 
vided by coimt rates for the wavelength range between a 
given lower limit and an upper limit -0.25 A. This cal- 
culation is performed by a summation of the energy fluxes 
times a weighting factor, WT, which has values of unity 
except for the first element, which has a value of 0.5. 

The weighting function is used to perform the energy in- 
tegration as in Figure C-2. The- value of X is set as a 
vector of 16 zeros by step 40. 

Steps 41 - 44 add -the value of X to GV, element by element. The first 
time through this addition* does nothing , but subsequent 
times it serves to fill out each increment since dE/dX 
for the end points was taken to be 0.5 • dE/dX for the 
interior points. Step 42 defines X equal to a 16-element 
vector of. 0,5 • dE/dX divided by the count rate vector 
and adds X to GV, element by element. Then GV is lin- 
early averalged. Step 43 is equivalent to 

I GV, 

AVE.ijg- 

The standard deviation, SIG, is determined by step 44, 
which is equivalent to 


Steps 33 - 
Steps 36 - 


SIG 




1/2 
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YP/3P 



UPPER LIMIT 

WAVELENGTH 


Figure C-2. Schematic Illustration of the Integration Technique 
Used for Calculating Energy Intervals 




Steps 45 - 49 I£ the value o£ FLAG = +1, the program branches to step 

'72. I£ FLAG = 0, the program continues to step 46. Step 
46 adds the new element SIG to the old vector SMAT, and 
step 47 adds the new element AVE, to the old vector AMAT. 

• Steps 48 and 49 add the new elements J and I to old vec- 
.tors XJ and XI. 

Steps 50 - 52 test the values of J and I (the upper and lower limits) , 
incrementing and looping if necessaiy. Step 50 loops to 
step 41 until all upper limits are calculated. Step 51 
resets J becaiise the lower limit is now changed. Step 
52 loops to step 36 until all lower limits are calculated. 

Steps 53 - 57 determine the minimum relative error (SIGM) and the cor- 
responding wavelength indices I and J. Step 53 divides . 
SMAT (a 25-element vector) by AMAT (another 25-element 
vector) , element by element. Step 54 (right to left) 
determines the minimum value of X and calls it SIGM. 

.Step 55 determines the index of X corresponding to SIGM, 
and steps 56 and 57 define the corresponding components 
of XI and XJ as the indices .sought. 

Steps 58 - 71 test the values of I and J and take appropriate action. 

Step 58 shunts the program to step 63 if the current 
value of SIGM equals that determined two iterations ago 
(stored in SIGM) or if I = 10 and J = JO . Otherwise , the 
program goes to step 59 . Step 59 branches the program to 
step 69 if SIGM equals the SIGM of the previous iteration 
or I and J do not occur on edges. Otherwise, the program 
goes to step 60. Steps 60 and 61 jump 10 and JO to values 
equal to the initial values plus twice the difference be- 
tween I and 10 or J and JO. Step 62 branches the program 
to step 30. Step 63' causes a printout of a 5-by’-5 matrix 
of relative errors. Step 64 prints a title for the matrix. 
Step 65 causes the terminal to skip a line. Step 66 prints 
out the optimum passband. In cases where one limit lies 
near 0 or 20.25 A, the optimal limits are not the same as 
the ones mentioned in the matrix title. Step 67 sets the 
value of FLAG =1, and step 68 branches the program to 
step 36. Steps 69 - 71 set 10 = I and JO « J and branch 
the program to step 30. 

Steps 72 - 77 print out the results on the terminal. Step 72 prints 

out the values of the linearly averaged conversion factors 
and the corresponding standard deviation. Steps 73 and 74 
calculate the logarithmically ‘averaged conversion factor 
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and the corresponding multiplicative factor, and step 75 
prints them out. Step- 76 is a header; step 77 prints out 
the conversion factors for free-free and blackbody spectra 
and for power law spectra. 

Steps 78-81 provide a notation that the program is finished (step 79) 
and a computation of the CPU time used (steps 80 and 81) . 
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V TASSBAND 
Cl] TSTRT^^rXll 

C2] AVE*-CR^GV^SIG-^SIGM^1CY^\0 
C3] SIGMA-^^^FLAGt-Q 
C4] I7-«-0.2 55<i80 

[5] CENTER GM TUBE:' 

C6] GM^U 

C7] 'ENTER GM TUBE AREA, IN SQ, CM, 

C8] GM^-LO .5 + ClOOOOOOOOO )xGA/^^:r-(-r/CM-^G'i'fxf/x(l^ 

1.9 8625"8)xjf-(-D 

C 9 ] ' ENTER ESTIMATED LIMITS , LONER ( S . 2 54 ) , UPPER ( <2 04 ) ; ’ 

[10] Z-t-D 

[ 11 ] lO-^WxXllI 

[12] J0^W\XL2l 

[13] 1 2 5 10 20 50 100 

[14] 'TEMP, RANGE 1 -1 0 0//ILL . » OK?iY~IES,OR ENTER 7 COMP, V 
ECTOR ) : ' 

ri5'i Y^a 

[16] y<*(143.8vP/)o .ij 

[17] !<■*( ( (J<170 )xj) + (y>170)xl70) 

[18] FF-(-LO .5 + (l 0 0000000 0 )xFF^(80pl )o . xZ-ef /FF-f ( ( {^/* 

“2)o .x7pl)4Z 

[19] FF<-L0.5 + (10 0000 000 0)x5Fr(8 0pl)o . xZ-c-f ( f;?* 

*■ 5)0 .x7pl)^(J-l) 

[20] 'ENTER NAVELENGTHS OF LINEiS) (&. 254 » 5204 ): » 

[21] 

[22] 'ENTER LINE/ CONTINUUM RATIONS):' 

[23] Z-f,Z^D 

[24] J^l 

[25] LINESiFFlYiJl ; ]-e-FF[ J[ J] ; ] xZ[J] +1 

[26] B5Cr[er];]^B5[Z[er];]xZ[dT]+l 

[27] ->((pY)SJ-f-j;+l )/LJFF5 

[28] CR^iXYxlE 9)x((GZ^+.xFF) ,(G//+*xBB)»(GM+.x6) 2 

80 p(80pl),f/* 2)) 

[29] GM<-\0 

[30] SETiX^JO-IO 

[31] JB-^((X0<3)x 3) + (( J0>2)xjG)-(I-G>8)x(4-Z)x(Z<4) 

[32] J'G-e((e/'0>78 )x78) + ((JG<79)xe7G) + (J£?58)x(4-Z)x(Z< 

4 ) 

[33] 1^10-2 

[34] JWO-2 

[35] AMAT-*-SMAT-*-XI^Xj^\0 

[36] F0f/:Zy'f-"i+r+l(e7‘-r) 

[37] Y^qXY 

[38] P/r-^O.SjCZ-l )pl 

[39] GF-f-CCr/y+.xFFCZZ;]) ,(r72'+.xBB[ZJ;]),(l'/2’+.xi?(2.Z)p(Zpl ), 

f7[ZJ]*“2))^(7F 

[40] Z*<-16p0 
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[41] COLiGY^GY^’X 

[42] ^?7<-GF+Z-«-0.5x(^>^'[J*],BB[e7’i],l,I/[J-]*“2)4-C'i? 

[43] i4Fff-t-+/(77^i6 

[44] i?JG!-e-(+/( (<?7“47ff)*2)^15)*0 ,5 

[45] ->-FLkG/Em 
[462 SUkT-^SMA'S ,SIG 

[47] AUk T^-AMk T ^k YE 

[48] XJ^XJyJ 

[49] XI^XI,I 

[50] HiJ0+2)^J<-J+l)/C0L 

[51] J*-J0-2 

[5 2] H-(( jO+2)>J-<-J+l)/i?OP/ 

[53] X^iSMAfiAMAT) 

[54] SIGMA <-SIGMA , SIGM-^- l/X 
[5 5] X<rXxSIGM 

[56] I^XIIXI 

[57] J^XJIXJ 

[ 58 ] ->( ( ( I- 10 ) A ( J =J0 ) ) V ( SIGM=SIGMA [ “ 2+pSJGMA ] ) ) / OUT 

[59] ■^({(JptIO-2)AiI^IO+2)AiJ*JO-i-2)f^(.J^JO-2))viSIGM=SIGMAL 
“l+pSIGMAl ) )/3TEP 

[60] IO^IO+iI-IO)x2 

[61] J0W0+(J-J0)x2 
[6 2] -^SET 

[6 3] OUT: 5 5 piSMAT^AMAT) 

[64] ' RELATIVE ERROR ABOUT ' ; f/[^J[l 3]] ; ' AND ' 

;J^[JJ-[13]] i M.’ 

[65] »• 

[66] 'OPTIMUM PASSBAND: »;P/[J];» TO »;f/[J'];* ANGSTROMS 

[67] FLAG^l 
[6 8] -J-BOP7 
[69] STEP:IO-^I 

[7 0] JO-^J 

[71] -^SET 

[7 2] END: 'LIN -AVERAGE GAMMA: 'iAVE-,' TO A SIGMA OF 'iSIG;' 

I 

[73] AVE^*C(+/»GV)^16) 

[74] ,ffJ(?-f-*C(t/((®(?7^jl7S)*2))Tl5)*0.5 

[75] 'LOG-AVERAGE GAMMA: ';AVEi' TO A MULT, OF 'iSIG;'-.' 

[76] 'GAMMAS CALCULATED:' 

[77] 2 7 pGV; 1 2 pffy[15] ,G7[16] 

[78] ” 

[79] 'END' 

[8 0] ^-*-121 

[81] » CPU TIME : ’ ; {X-TSTRT ) ^- 6 0 ; » SECONDS , ' 
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APPENDIX D. CALCULATION OF CONVERSION FACTORS 
ASSUMING A KNOWN PASSBAND 


The purpose o£ the program called KPASSBAND is to determine the 
conversion factor for a GM tube or an ion chamber or a set of GM tubes 
or ion chambers when the passband is known or assumed. The program is 
outlined in Figure D-1 and listed at the end of this appendix. 

The operation of this program is as follows. The GM tubes are 
entered serially (i.e., as one long vector with one tube appended on 
the next) . The limits are entered, and the spectra are calculated. 
(Temperatures may be specified as well as lines.) The incident energy 
flux within the specified passband is calculated. 

The first tube is treated by calculating the count rate .for each 
of the 16 incident spectra. The corresponding conversion factors are 
calculated, and statistics are 3nm on these conversion factors. The 
next tube receives the same treatment, and so on. After the last tube 
is treated, the program prints out the average relative error (linearly 
and logarithmically) and the resulting conversion factors and uncertain- 
ties. Each conversion factor may be printed out, if desired. 

In detail, the program operates as follows. 

Steps 1 § 2 set initial variables. 

Steps 3 § 4 call for arid enter the GM tube efficiencies . 

Step S determines the number of GM tubes inserted in 

step 4. 

Steps 6 § 7 call for the entry of the GM tube area (assumed to 

be the same for all tubes), calculate 1/hc times the 
efficiency times the area, normalize it to an integer 
form accurate to one half part in 10^°, and define the 
normalization constant as NORM. 

Steps 8 § 9 call for the entry of the limits, in angstroms. 

Step 10 defines the vector, W, of wavelengths, in angstroms. 

Steps 11 § 12 determine the indices of W corresponding to the limits 
entered in step 9 . 
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STEPS 1-16 


STEPS 17-21 


STEPS 22-29 


STEPS 30-32 


STEPS 33-34 


STEPS 35-38 


STEP 39 


STEPS 40-44 


STEPS 45-47 


STEPS 48-52 


STEPS 53-56 


Figure D-1. Flowchart for F^PASSBAND Program 
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Step 13 determines a vector of the indices ranging from the 

lower limit, I, through the upper limit, J. 

Steps 14 - 16 pre “define Y equal to a range of temperatures and ask 
for the entry of the temperatures wanted. 

Step 17 defines Y = hc/WkT. 

Step. 18 limits Y to values <. 170 and then takes the exponent to 

the base e of Y. 


Steps 19 § 20 define the free-free and blackbody continuum spectra. 

These spectra are normalized to the largest value for 
each temperature and then converted to an integer form 
accurate to one half part in 10® . 

Step 21 defines the power law spectra in integer form. 

Steps 22 - 25 ask for the entry of the line wavelengths and cor- 
responding line/ continuum ratios . 

Steps 26 - 29 form a loop which replaces the old spectra values 

(free-free and blackbody only] at wavelengths which 
include lines with the old spectral flux value times 
the sum of the line/continuum ratio and 1. 


Steps 30 - 31 define a weighting function, WT, such that the two 

end increments of the energy flux integral are weighted 
one half the weight of the other intervals, and define 
an energy flux vector (16 components, one for each spec- 
tirum) using the spectra and the weighting function. 

The expressions within the parentheses are equivalent 
to (using the leftmost expression as an example); 

(WT +,x FF[XY;]) = I WTj ♦ FF^j 

Steps 32 - 39 form a loop which calculates the statistics for each 
individual tube. Step 33 defines GM equal to the 80- 
component efficiency vector of the Jth tube, times the 
wavelength, and converted out of normalized units. 

Step 34 defines the vector GV equal to the energy flux, 
EFLUX, divided by the count rate (the expressions in 
parentheses) and appends this vector onto the vector 
, GAM. Steps 35 and 36 calculate the logarithmic average 
and standard deviation, and steps. 37 and 38 calculate 
the linear average and standard deviation of the resul- 
tant conversion factors. Step 39 increments J, tests it, 
and branches the program to step 33 if J is less than or 
equal to the number of tubes . 
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Steps 40 - 44 print out the results o£ the above calculations . 

Steps 45 - 52 inquire i£ all calculated conversion £actors are 
wanted and, if not, branch to step 53. Otherwise, 
the program prints out all conversion £actors cal- 
culated for each tube. 

Steps 53 - 56 announce the end o£ the program and print out CPU time. 
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V KFASSBAND 
[1] TSTRT-^Xll 

121 GM^GAM^SIG-^-LSIG^A VE^LA VE^EFLUX^^rGV^XX ^ 1 0 

C 3 3 * ENTER GM TUBES ; ' 

[43 GUTB-^U 
[53 N-<ripGMTB)iBO 

[63 CENTER GM TUBE AREA, IN SQ. CM. 

[73 GMTB^l 0.5 + ( 10000000000 )x GMTBiNORM*- f / GMTB^GMTB x ( 1* 

1 .986 2S’”8 )xZ^O 

[83 CENTER LIMITS, LOWER .25A) ,UPPER(.^20A) i ' 

[9.3 M 

[10 3 25x180 

[113 J-«-f/iJ[l3 

[123 J’-6-I7iJ[23 

[133 XX^~lfI+\il+J~I} 

[143 If’ 1 2 5 10 20 50 100 

[153 'TEMP. RANGE 1-100 MILL. » OK? (I=Iff5, OR ENTER 7 COM 
P. VECTOR)'.^ 

[16 3 If-D 

[173 Jf'(143.84-J/)o,iy 

[183 If-*(( (Ial70)xl70) + (I<170)xI) 

[193 Fff-L0.5 + (100000000)xPF^(80pl)o .x^-.-r/J’F-^CCP/* 2)»,x 
7pl)^'I 

[20 3 SS-f-LO .5 + (10000 0000 )xBF4(80pl . xl-t-f /SBf-( (F/i^S)" . x 
7pl)v(J-l) 

[213 Zf-Js? 2 80 p(80pl ) ,IvIf-L/If-f/* 2 
[2 23' 'ENTER LINE WAVELENGTH (.S') (A)i' 

[23 3 If-.I^F/iIf-D 

[243 ’ ENTER LINE / CONTINUUM RA TIO ( B ) ; ' 

[25] If.If-D 
[2 6 3 «7’f-l 

[27] LINES t FFL lUl ; ] f-FF[ I[«73 ; 3 xl[*^3 +1 
[283 BB[l[«r3;3f-B5[l[«/’];3xl[J3+i 

[293 *((0Y)^J^J+1)/LINE3 

[30] J/ff-O.S .((~2+pII)pl) ,0.5 

[31] SFB(/Zf-(f/2'+.xFF[II;3).(r/2’+,xBB[yj;3),((72’t.xZ[II;3) 

[323 J'fl 

[33] (J-l)x80) + \803xI/xi70i?Mx 

1B"10 

[343 GAM^GAM , GV^^rEFLUXi ( ( . xFF ) , ( GM+ . xBB ) , ( GM-i - . x Z ) ) 

[353. LAVE*-LAVE,*(i-/ (^GV)^ie) . 

[363 LSIG-<r’LSIG,*((+/ ((a(GV^LAVEtpLAVEl))*2)^A5)* 

0.5 ) 

[373 AVE-^AVE,(+/GVrlS) 

[38] SIG^SIG .(+/(( GV-A VEl p4 7B3.) *2)415)* 

0,5 

[39] -*(N>J-^H^1)/ RERUN 

[403 'OVERALL UNCERTAINTY i LINEAR RELATIVE ERROR-' i-i-/ (SIGi 
AVE)^N',' LOGARITHMIC SIGMA -';+/( LSIG^N) ; * . ’ 
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[41] 'RESULTIIIG COEVERSIOE FACTORS ARD UNCERTAINTIES:^ 
E421 ^({5jJ)p(xN) ,AVE,SIG^LAVE,LSIG) 

!;43] • LIN GAMMA., StGMA LOG GAMMA, SIGMA^ 

[443 ’• 

[453 ' PRINT OUT ALL GAMMAS ? ( 1 ^YES ) ; ’ 

[461 r-e-0 

[473 /END 

[483 

[493 NXTi'TUBE ' ;J ; ' : ' 

[503 2 7 pC?/lW[ ( (c7’-l)xl6) + ll43 ; 1 2 pGAMC((V-l)xl6)i- 

1 53, GAME 6 2 

[513 '» 

[521 -yCN^JW+l)/NXT 
[533 ENU:' ’ 

[541 ’END^ 

[5 53 ;f-f-X21 

[563 ' CPU TIME : ' ; ( X-TSTRT )-^ 6 0-; ’ SECONDS , ' 
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APPENDIX E. RESPONSE BEHIND A CONICAL APERTURE 


This appendix describes a program designed to determine the 
relative response o£ a GM tube' (or ion chamber) located behind a 
conical aperture. The program includes calculations determining 
the shadowing of the GM tube window by the aperture. The name of 
the program is CONE. The program called CONE IN generates the input 
for CONE. 

The CONEIN program is flowcharted in Figure E-IA. CONE is flown- 
charted in Figure E-IB. The GM tube absorption and efficiency are en- 
tered with the GM tube geometry (window area and gas path length) , 
spectra are determined for future use, the aperture geometry (Figure 
E-2) is entered, and the needed angles are calculated. The overall 
approach involves placing a 6-by-6 grid of squares over the face of 
the GM tube. The squares are dimensioned such that the four comer 
squares lie outside the radius of the window (Figures E-3 and E-4) . 

The area and the center of gravity for each square are calculated. 

The aperture -shadowing problem then becomes a problem of determining 
how much of each areal element is left unshadowed and what the approxi- 
mate center of gravity is, so that the efficiency of each element can 
be calculated, weighted by the area of that element, and then all ele- 
ments summed to determine the final effective efficiency. 

The program CONEIN is described in detail as follows . 

Step 1 clears a number of variables . 

Steps 2-6 require the entry of the gas absorption of the GM 

tube and the efficiency of the GM tube and calculate 
the window transmission. 

Steps 7-11 determine values hc/XkT. Step 7 sets up the standard 
temperature range required for step 9 . Step 10 deter- 
mines hc/XkT, and step 11 checks the maximum limits 
and then exponentiates. 

Steps 12-24 require the entry of the GM tube geometry and the 
aperture geometry (Figure E-2) and define variables 
R1 and R2. 

Steps 25 - 26 require the entry of the off-axis angles wanted. 

Step 27 determines the off-axis angles needed (Figure E-5) . 

Steps 28 - 30 determine the -areal elements and the dimensions shown 
in Figure E-4. 
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STEPS 1-6 


STEPS 7-n 


STEPS 12-16 


STEPS 17-26 


STEP 27 


STEPS 28-41 


Figure E-IA. Flowchart for CONEIN Program 
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STEPS 1-15 
STEP 1^ 

STEPS 17-78 
STEP 79 
STEPS 80 & 81 

STEPS 82-85 

STEPS 93-97 

STEPS 98-100 
STEP 101 

STEPS 102-104 
STEPS 105-118 

STEP 1 19 

STEPS 120-125 


Figure E-IB. Flowchart for CONE Program 
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Figure E-4. Dimensions of Grid Network Used in CONEIN and CONE 
Programs 
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VOLUME (CYLINDRICAL) 


Figure E~5. Side View Illustrating Aperture Geometry and 
Defining Angles 
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Steps 31 - 36 deteimine the X and Y (Figure E-4) coordinates of the 
center of gravity of the areal elements . 

Steps 37 - 39 determine the differentials of the X and Y center-of- 

gravity matrices (i.e., the distance betvreen the center 
of gravity and the edge of the areal element) . - 

Steps 40 - 41 serve to define variables and announce the end of the 
program. 

The operation of program CONE is as follows. The window of the 
GM tube is overlaid with a 6-by-6 grid of squares such that the four 
comer squares are outside the window proper (Figure E-4) . The amount 
of shadowing of each square is individually determined, and the shift 
in the -center of gravity of the area of each square is approximated. 

This effective center of gravity is later used to estimate the effec- 
tive gas path length. Because of the symmetry of the conical aperture 
and cylindrical GM tube system, only half of the areas of the squares 
need be determined (Figure E-4) . The uniformity of the incident flux 
of X-ray photons across the face of the aperture is tacitly assumed. 
Thus, the product of the area and the efficiency is actually calculated. 

In detail, the program CONE operates as follows . 

Steps 1, 2, 121, 122, 123, 124, and 125 generate statistics on central 
processor time (typically 100 seconds on an IBM 360/50) , 
turnaround time (typically the order of minutes) , and 
the number of users on the system. 

Steps 3-9 set initial values for variables. Step 5 determines the 
value of the area times efficiency, AXE, were the GM tube 
viewing the source head-on with no aperture. Step 6 de- 
fines the wavelengths used. Steps 7 through 9 refer to 
geometric constants (defined in Figure E-5) . 

Steps 10 - 15 list the aperture dimensions used for this particular 
run and ask if stopping the program or obtaining de- 
tailed output is desired. The geometry and the dimen- 
sions are defined in Figure E-5. 

Step 16 routes the program to step 98 if J = 0. This is done 

to determine the relative count rate (a 16-element vec- 
tor, one component for each of 16 spectra — seven free- 
free, seven blackbody, and two power law) if the GM tube 
is viewing the source head-on. 

Steps 17 - 19 allow the addition of input to the program if additions, 
were requested. 



Steps 20-22 determine offset values (Figure E-5) . 

Step 23 sets the relative area matrix, AR, equal to a matrix 
of zeros . This sets up the program for the condition 
of total- shadowing. 

Step 24 determines if the aperture is totally shadowed. This 

condition occurs when the projection of either aperture 
on the plane of the GM tube window does not overlap the 
window. (Figure E-6) or when the projections of both 
apertures on the plane of the tube window are mutually 
exclusive (Figure E-7) . If the window is completely 
shadowed, "Ae program branches to step 79. 

Step 25 redefines the relative area in terms of no shadowing. 

Step 26 determines if the window is shadowed (Figure E-8j . If 

the window is not shadowed, the program branches to step 
79. 

Steps 27 - 29 set values for variables. 

Steps 30 § 31 determine if only the lower mask, Rl, effectively masks 
the tube window (Figure E-9) ; if that condition is true, 
step 31 branches the program to step 38. 

Step 32 branches the program to step 47 if only the upper mask 
effectively shadows the tube window (Figure E-10) . 

Steps 33 - 36 determine if the intersection of the projections of the 
two masks on the plane of the window lies outside the 
radius of the window (Figure E-11) . Steps 33 and 34 
calculate XI and X2 from the law of cosines (Figure E-11) . 
Steps 35 and 36 test whether XI ;< X2, and, if so', branch 
the program to step 47. This series of steps differs 
from step 31 in that it tests for different conditions. 

Steps 37 § 38 calculate the angles XI and MAX, defined in Figure E-12, 
using the law of cosines. 

Step 39 is similar to step 36. It branches the program to step 

47 when XI ^ MAX; i.e., the mask of radius R2 does all 
of the effective shadowing. 
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IF (OF2+R)<R2 AND 
(OFl+R) <R1 
GM TUBE WINDOW IS 
UNOBSTRUCTED. 



Figure E-8. Condition for No Shadowing 



Figure E-9. Condition for Only Lower Mask to 
Shadow the GM Tube Window 
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Figure E-12. Geometry Defining Angles XI and MAX 
(Steps 37 and 38) 
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Step 40 determines the angles between the offset line, OFl, and 

radii dravm to the intersections of the grid lines .and 
the edge of the tube window (Figure E-13) . The expres- 
sion to- the right of the comma deteamiines the angles of 
the intersections with the horizontal grid lines; the 
• expression to the left of the comma deteimiines the angles 
of the intersections with the vertical lines. 


Steps 41 - 43 order the angles in ascending order, truncate the vector 
of ordered angles so that it includes only angles less 
than MAX, and then truncate angles less than XI. 


Steps 

44 

Step 

46 

Step 

47 

Step 

48 

Step 

49 

Step 

50 

Step 

51 

Steps 

52 

Step 

54 

Step 

55 


§45 define vectors with the number of components equal to 
that of the vector ANG and of values of the radius R1 
(RAD) and OFl (OFST) . 

branches the program to step 58 if the mask of radius 
R1 is the only mask shadowing the window. 

is the same as step 46 and covers the geometry shown 
in Figure E-14. 

is the same as step 40, but the upper mask (of radius 
R2) is used rather than the lower mask (Figure E-13) . 

orders the elements of vector A2 in ascending order. 

calculates the angle MAX illustrated in Figure E-15. 
This is the angle at which the outer edge of the window 
intersects the projection of the upper mask. 

truncates the vector A2 so it includes values less than 
and equal to MAX. 

§ 53 branch the program to step 55 if only the upper mask 
shadows the window. 

truncates the vector A2 to values less than or equal . 
to X2. 

redefines the vector RAD. The old vector RAD contains 
components of values equal to Rl, with each component 
corresponding to an element in the vector ANG. RAD is 
appended on the end of a vector of components equal to 
R2, each corresponding to a component of A2. The re- 
sulting vector is called RAD. 
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FIRST- TERM SECOND TERM 



Figure E-13. Geomeiry Showing the Intersections Between 
the Grid Lines and the Profectior) of the 
Lower Mask 
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Figure E-14. 


Geometry for Shadowing by Only Lower Mask 



Figure E-15. Definition of the Angle MAX as the Inter 
section of the Proiection of the Lower 
Mask and the GM Tube Window 
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Step 


56 


Step 57 
Step 58 


Steps 59 
Step 61 
Step 62 

Step 63 
Step 64 
Steps 65 

Step 67 

Step 68 
Steps 69 

Step 72 


operates similarly to step 55 except that corresponding 
values of OFl or 0F2 are used. 

determines a new vector of angles , ANG, consisting of 
old vector ANG appended to the end of the vector A2. 

calculates a vector of values consisting of segments of 
circles (Figure E-16) . These segments correspond to the 
angles and radii contained in ANG and RAD, respectively. 
The equation for the area of a segment, seg, is 

seg = 1/2 • (r^) • ((2 • half-angle) 

- sin (2 ♦ half-angle)) 

8 60 calculate the coordinates of the intersections of the 
grid and radii R1 and R2 based on the system shown in 
Figure E-4. 

branches the program to step 64 if the largest value of 
Y > (D - B) . This means that if only grid 1,1 was shad- 
owed (Figure E-4), the. program does not branch to step 64. 

sets the area of grid 1,1 equal to the totally unshadowed 
value, ARO (1;1) , minus the value of segment (2) . This 
result follows because seg (1) =0, and seg (2) corre- 
sponds to the segment that crosses the edge of the window. 

branches the program to step 79. 

sets the index K = +1. 

8 66 define lY equal to the value of Y(l) divided by D, rounded 
off to the next higher integer, and branch the program to 
step 68 if lY = 1. 

sets the values of all components of the area matrix, AR, 
which are completely shadowed equal to 0. This step re- 
quires that lY be greater than 1. 

branches to step 75 if ANG(K) = X2. 

- 71 define lY and IX as Y(K + 1) and X(K + 1) divided by D and 
rounded off to the next higher integer and, if IX = 3, 
branch the program to step 73. 

sets the areas of grids lying outside the grid in ques- 
tion, but on the same row, equal to 0. 
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Step 73 


defines the dummy variable MAX as the area of a grid 
(times 2, since by syinnietry only half the grids need 
be calculated) between its maximum Y values (the les- 
ser of either the outer limit of the window (lY • D or 
R + (3 • D)) and the value of Y(K + 1) provided that 
the present value of K + 1 does not indicate the tran- 
sition case in which the shadows from the two ma^ks 
cross and provided that K + 1 is not the last angle. 
Otherwise, MAX = 0, 

Step 74 determines the area of the grid. This area equals the 

difference between the areas of the segments K and K + 1 
(resulting in a strip across the face of the tube) run- 
ning perpendicular to the Y axis) , ' plus the prior value 
of the grid if ANG(K) ^ X2, plus MAX, and minus a rec- 
tangle laid on the strip which has a height equal to 
(Y(K + 1) - Y(K)) and a length equal to (IX - 1) • D ♦ 2. 
This step leaves, the areas shown in Figures E-17 and E-18. 

Step 75 increments K and, if K is still less than the number of 

components in the vector ANG, branches the program to 
step 68- Advancing to the last intersection drops the 
program out of this loop. 

Step 76 defines MAX as a vector of two angles (Figure E-19) . 

Step 77 sets MAX as a vector of the corresponding areas of the 

segments of the previous vector MAX. 

Step 78 determines the area of the last grid which needs calcu- 

lation. This is the grid which contains the intersection 
between the shadow and the edge of the window'. This area 
is equal to the previously calculated area, plus the dif- 
ference between the two segments calculated as MAX (abso- 
lute value) , minus the area of rectangle contained in the 
difference of the two segments. See Figure E-19. 

Step 79 determines the variable L to be a 6-by-5 matrix of I’s, 

Steps 80 § 81 branch the program to step 93 if J = 0 or TX(J) = 0. 

Steps 82 8 83 determine approximate center-of-gravity matrices for the 
Y(CGY) and X(CGX) directions using initial centers of 
gravity (CGY, CGX) and a fraction of the difference 
between the center of gravity and the edge of a grid 
(CGYD, CGXD) . A simple linear change was used as an 
approximation in Figure E-20. 
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Figure E-19. Definitions of AAAX 1 and MAX 2 (2 Cases) 
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Figure E-20. DefiniHons of Points of Centers of Gravity 
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Step 84 defines L as a 6 -by -3 matrix of the effective gas path 

lengths, in units of LO, (See Figure E-21.) 

Step 85 sets L equal to the minimum of either L defined in step 

84 or the secant of TX(J) . (See Figure E-22.) 

Steps 86-91 determine that if optional program stopping is desired 
- the program stops (steps 86, 87, 88), asks what is de- 
sired, and then goes on. If stopping is not desired, 
the program branches to step 92. 

Step 92 erases variables no longer needed. 

Step 93 defines X as a vector (each element corresponding to a 
different wavelength) of the window transmission. The 
window transmission is defined as the initial transmis- 
sion to the secant TX(J) power because the change in 
path length occurs in the exponential argument of the 
expression for the transmission (i.e., T = exp (p/p)px 
sec TX(J)) 

Steps 94 - 96 define the effective area X efficiency, AXE. This is 
equal to the sum, for each grid, of the window trans- 
mission (X, defined in step 93), times the area of the 
grid, times 1 minus the gas transmission to the L power. 

Step 97 corrects for area foreshortening due to off-axis look- 

angles . 

Steps 98 - 100 calculate a 16-component vector of the' relative counting 
rates for free-free, blackbody, and power law spectra. 
Step 98 is equal to 

CR = Z AXEj • Wj • Tjj 

where 
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X 



y 

Figure E-21. Geometry for the Effective Gas Path 
(One of Two Cases) 
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Step 101 
Step 102 
Step 103 
Step 104 
Step 105 


branches the program to step 104 i£ J 0. 

defines CRO as equal to CR , ' 

branches the program to step 119 . 

defines the relative counting rate, RC, as CR/CRO. 

sets the Jth component of the vector RCR equal to the 
average value of RC, or 

I RCj 

I 16 


Steps 106 - 118 branch the program to st^ 119 if the option (steps 14 
and 15) is not desired; otherwise, the program asks if a 
detailed output is desired. If the answer is 1 (yes), 
the program prints out a list including the angle TX(J) 

(in degrees) , the relative area (coijpared to looking head- 
on), the relative counting rates determined for each of 
the 16 spectra, and the linearly and logarithmically aver- 
aged relative count rates and their standard deviations . 

Step 119 branches the program to step 17 if, after incrementing J, 

J is not larger than the number of angles. 

Steps 120 - 125 announce that the program is finished and determine time 
statistics . 
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V CON BIN PUT 

[1 ] RCR-*-TX->rAEOO-*-CGXO-^CGyO-irCGXD-<r-CGYD<-\0 

[23 'ENTER GAS ABSORPTION:' 

[33 TRGAS-i-l-X-i-a 
[43 'ENTER GM TUBE:' 

[53 X^[l 

[63 TRWIN-^X ^1~TRGAS 

[73 J-6- 1 2 5 10 20 50 100 

[83 'TEMP, RANGE 1-100 MILL, « K OK? (X=XESyOR ENTER 7 CO 
MPONENT VECTOR ) : * 

[93 

[10] 2’-<-(143.8f(0.25xi80))o 

[11] 2’'e-*(( (2’<170)x2’) + (2'>l70)xl7 0) 

[123 'ENTER AREA OF GM TUBE iSQ , CM,):' 

[13] 

[14] J?^(ZtO1)*0.5 

[15] 'ENTER LENGTH OF GM TUBE (CM.);’ 

[16] LO-^U 

[17] 'ENTER RADII OF TWO MASKS iCM,), ASCENDING ORDER:' 

[18] Z-(-D 

[19] i?l-t-Z[l] 

[ 20 ] R2^Xl2l 

[21] 'ENTER CORRESPONDING HEIGHTS OF MASKS iCM,):' 

[22] Z^D 

[2 3] HTSl*-XZll 

[24] HTS2-<rXL2'} 

[25] 'ENTER MAX, THETA, GRID (DEGREES , DEGREES) : ' 

[26] Z-«-(O^180)xi:-(-a 

[27] TZ-«-0.(Z[23xi(LZ[1]tZ[2])) 

[28] AROO-^ 6 3 p.v8 

[29] AROOL \ 3 ; 33-t-4>Ai?CO [ 1 ; l^AROOL 6 ; l^AROOl 3 + t 3 ; 33+- 
0.096043 0.046656 0 

[303 i4i?CC-{-L0. 5 + (10 000 00)x^i?CO 

[313 X*- 0.48032 1.3577 0 0.5 1.5 2.23265 

0.5 1.5 2.3861 
[323 CGXO^ 6 3 pZ 
[333 CGZC[3 + i3;]-*-Jsj(})isjCffZC[ \3;] 

[343 CCIC-*- 6 3 p( ,4>6}(CCZC[ i3;3 + 3)) 

[351 CGZC[ i3;]-<-6-!s?<j)S?CCZCC3 + i3;] 

[363 CCIC[ 1;3>CCJC[ 6 ; 33-6-0 

[37] CGXD^CGX0-(6 3 p(0 1 2)) 

[38] CCZ2?[1 ;3]'e-CGZZ?[6 ; 33-6-0 

[39] CffJP-6- 6 3 p 0.3861 0.2327 0 0.5 
0.5 0.3577 0.5 0.5 0.48 0.5 0.5 
0.52 0.5 0.5 0.379 0.401 0.249 

0 

[40] RCR-^QxTX 

[41] 'SAVE THIS WORKSPACE, LOAD "CONE", COPI "INPUT INC 
ONE' ' , AND RUN ' ' CONE ' ' . ’ 

V 
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V CONE 

Cl] ■ TTSTET^X20 
C2] TSTRT^X21 
C3] J<-0 

C 4 ] A NG-^X^I^EA D^OFST^xO 

C5] AXE*-iOR*2)xTEi/INx.l-TRGAS 

[6] r/-f-0.2 5xi80 

[7] 2?'H(2*~1.5)xi? 

[8] i4-«'0.2 92 89xi? 

C9] B-^-0.2203xi? 

[ 10 ] 'ffW TUBE OF RADIUS CM., LENGTH CM., AND 

AREA SQ.CM. ^ 

[11] ^MASKS OF »;i?l,i?2;’ ClI . RADII AT HEIGHTS OF 'iHTSl, 
HTS2; ^ _CM . ’ 

[12] MAX<r{{ 3o( (J?+i?2)^7/2’5’2))[(“30((i?+i?l)f7/2>£'l)))x 
1804O1 

[13] 'MAX. THETA = * ;2'X[ pffZ] xi 80 ^ Ol ; ' DEG.\MAX VISIBLE = ’ 
iMAXi' DEG.' 

[14] ’ OPTIONS • UA NTED ? ( 1 =IES ) ’ 

[15] OPTN*-U 

[16] -^U = 0)/CRCALC 

[17] CY CLEA NGi-^i OPTII^ 1) / CY CLEA NG ^ 3 

[18] ’ i' ^OPTION NOUi' 

[19] X^U 

[ 20 ] OFI^HTSI^ZOTXIIJ 

[21] OFl^HTSl^ZoTXW^ 

[2 2] OFZ^^OPl-OFl 

[2 3] 4i7-e- 6 3 pO 

[24] ^( ( ( OFl -R1 )^i?)v ( ( 0F2-JF2 )>i?)v ( (i?l+J?2)S0F3) )/ LENGTH 

[ 25 ] AR-^ARO^(oR*2 )xAROO^(+/ + /AROO) 

[26] ->( ((0F2+R)^R2)x( (OFl+R) SRI) )/ LENGTH 

[27] A2^A NG^RA D-<-OFST^ \ 0 

[28] ;s72-^00.5 

[29] .Yl-^0 

[30] i?1077Ly^((0F3+i?l)^i?2) 

[31] -^RlONLY/SKIPl 

[3 2 ] ■*{iqFl+R)SRl)/SKlP 2 

[3-3] Y2-^“20( ( (F2*2 ) + (OF3*2 ) -i?l *2 ) 4 2 xi?2xOF3 ) 

[34] Yl-«-"2o(( (F2*2) + ((?F2*2)-F*2)42xF2x£?F2) 

[35] R20NLY^X1<X2 

[36] ->R20 NLY/SKIF 2 

[37] Yl-<-‘*2o(( (F2 x2oY2)-(?F3)47?1 ) 

[38] JKIPl iMAX^~2o(^ ( (F1*2) + ((7F1*2 )-F*2)42x7?ixi7Fl) 

[3 9] -^{MAXSXD/SKIPZ 

[4 0] ^//G^(“1o((0,Z),2xZ?) 4F1) )."20(( (I?xiL(Fl-(FlOFl ))4D)t(i)| 

OFl ) ) 4F1 ) 

[41] ANG^ANGLUNGZ 

[42] A NG‘<r {{ANG <MAX) /ANG ) , MA X 

[43] ANG-^XlAiANG>Xl)/ANG) 
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[iJ-43 RAd-^ipANG)QRl 
[45 3 OFST-^ipANG)pOFl 
[463 ^RIORLI/SEGMNT 

[47] SKIR2i-*i{0F2-^R)^R2)/SEGMRT 

[48] A2-i-("loC(0,Z),2xZ))^i?2)) ,“2o(((PxiL(i?2-(Z)!0j?2))'fZ?) + (Z)| 
0F2))iR2) 

[49] 42-<*i42[4423 

[50] MAZ-f-“2o(( (i?2*2) + (OF2*2)-i?*2)T2xi?2xi?F2 ) 

[513 A2-^{{A2<MAX)/A2),MAX 

[52 3 -yR20NLX/SKIPZ 
[533 ■^iX2>MAX)/SKIPZ 
[543 A2~(-QiA2<X2)/A2) ,X2 

[553 SKIP3tRAD*-UQA2)pR2) ,RAD 

[56] OFST^i(.pA2)pOF2) »OFST 

[57] ANG*A2^ANG 

[58] SEGMRTiSEG-<-0 . 5 x ( RAD*2 ) x ( ( 2 x///(?) -lo( 2 xAF(? ) ) 

[59] Y^(,3^D)-fOFST-RAD>^2oANG 

[60] X^RADxloANG 

[61] HyipANGl>(.D-B))/KSET 

[62] Ai?Cl ;1]^/Ii?0[l;l3-5FG[2] 

[63] ^LENGTH 

[64] KSETiK^l 

[65] IX^-YlLlliD 
[663 Ml=IY)/CALC 

[67] >ii?[xCrr-l);]-«-((TY-l).3)pO 

[683 CALC:HANGLK2=X2)/TCT 

[69] IXHYlK-nl^D 

[70] JX't-rZ[if+l]4P 

[71] ->iIX=3)/3TCT 

[72] Aif[JX;JX+i (3-JX)3^(3-rX)p0 

[73] 55'Cr:M4Z-t-2xZ}x( (( (lJxZ})L(i?+3xi?) )-X[X+l3)x(4ii/G[i!:'+l];^Xl )) 

x( (K+l)^ipANG) ) 

[74] ARLIY ilXl^(SEGLKi-ll-SEGlK2 ) + (4i?[ JI ;TX] x(;l/7G[F] =X2 ) ) + 

M.4X-CJ[X+l3-X[is:3 )x2xZ?x( ) 

[753 TCT:HipArJG)>K^K-H)/CALC 

[763 /.MX^" 20 ((( |{Y[X3-3xZ))).( | ( { (IJ- 3 ) xp ) [i?) ) ) vi? ) 

[773 AM X-<- 0 . 5 X ( i? * 2 ) X ( ( 2 X /M X ) - 1 0 ( 2 X AM X ) ) 

[78] 4i?[JX;IX3-e-Ai?CJJ; JX] + ( 1 (/MX[13-AMX[2]))-2 xDx( (( JJxZJ)[( 
F+3 xD))-X[X3)x(JX-1) 

[79] LENGTH:L-<r 6 3 pi 

[80] ■*(J^0)/NEWEPF 
[813 ■*(.TXLJ1=0)/NEWEFF 

[823 CGY^Dy^CGY^CGYO+CGYDytil-iARURO)) 

[833 GGX-f-Dx GGX-J-GGXG- GGXZ)x ( 1 ->li? ) 

[84] L-f-( (i?x2o("lo(GGX^2?) ))-(GGJ-3xp) )^(£Oxlo2’X[<7] ) 

[853 L<-Llii-20TXLJ1) 

[863 -^iOPTN^D/NEP/EFF-l 

[87] ^PRINT iAR^ARO iL)?il=YES} ,ST0Pi-^^YES) , GO ON (9=JFS):» 
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[8 83 M 

[893 ■y(,K^l)/NE]lEFF-l 
[903 ^AR^ARO^ 

[913 »5';L 

[9 23 'M^ANG^SEG^RAD^CGX^CGY^xO 
[933 imJEFF : X^TRl/IN* ( 4 20ffJ[/3 ) 

[9 43 AXE^-i-ZdiX^ .^ARLilD^il-iTRGASo 

[9'53 AXE-^AXE+i + /iiXo :xARL%2'i)>^(.l-(.TP.GAS<^ ,*TA. ;23))) ) 

[963 AXE-^AXE-i-i+/ {{Xo .>^ARL i^l)>^il-(.TRGAS° .*LL ;33) )) ) 

[973 AXE^AXExlOTXZJl 

[983 CRCALCiCR^UXEi^.') + .x( v?) 

[09 3 Cif-eCi?, C(/l.Y27xj/*“4) + .x(f2’-l )) 

[10, OJ CR^CR,i(AXE^:/)-i-.x(^ 2 80 p ( 80p 1 ) ,Jy’*-2 ) ) 

[1013 ■^U^O'i/SKIPSET 
[1023 CRO-irCR 
[103 J -^RUNON 
LlomSKIPSETtRC^rCRiCRO 
[1053 ReRZJ2*-+/RC-iiQ 
[1063 HOFTN^D/RUNON 

[1073 ^DETAILED OVTPUT RAETEV? {l-YES) i ' 

[1083 Y*-D 

[1093 -y(Y^l)/RUN0N 

[1103 'AT THETA = ’ ;2’X[J’3 >'18 0^01 ; • DEGREES,' 

[1113 ' EEL, AREA t * (+/+/yli?)v(o/?*2 ) )x2o2’X[J3 

[1123 •^iX=Q)/RUN0n 

[1133 ’ REL. COUETIHG RATES i 2 7 pRCi 1 2 pCM 
153 ,i?C[163 ) 

[1143 X-t-+/RCTl& 

[1153 ’ LlE-AVEi +,-* ';(+/(( (i?<7-A:)*2 )vl5) )* 

0.5 

[1163 X^*( (+/©i?C’)^16) 

[1173 ' LOG-AVE: 'iXi' TO’ A FACTOR OF ' ;*(i+/ (ii®RC^X)* 
2)^15))*0.5) 

[1183 * ’ 

Lll9lRUN0NiHipTX) ^J-t^J+1 ) / CYCLE A EG 
[1203 'FIEISHED' 

[1213 ^-«-l21 

[1223 'CPU TIME iiX~TSTRT)i60i' SEC.' 

[1233 X^120 

[1243 'REAL TIMEi' i(.X-TTSTRT)*6.0;' SEC.' 

[1253 X^l23i' USERS.' 
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APPENDIX F. DETERMINING TEMPERATURES FROM COUNTING RATE RATIOS 


The program called RATIOS can be lased to determine the £ree-£ree 
and blackbody. temperatures corresponding to either count rate ratios 
or energy flux ratios providing that the energy flux ratios are derived 
by multiplying count rates by constant conversion factors. The program 
is flowcharted in Figure F-1 and listed at the end of this appendix. 

The operation of the program is as follows. The- product of the 
efficiency and the area of each GM tube is entered, multiplied auto- 
matically by X/hc, and converted to an integer foamat accurate to one 
half part in 10® . If energy flvjx ratios are to be used, the product 
of the efficiency (an 80-coraponent vector) , the area, and the conver- 
sion factor used is entered. The program calculates the ratios ob- 
tained for 18 temperatures ranging from 0.4 to 100 million °K. Due to 
core limitations, these calculations are run in blocks of six tempera- 
tures. The spectra for six temperatures are calculated, and then the 
ratios are calculated. The next block of ten^eratures is handled in 
the same way imtil all three blocks of temperatures have been processed. 
Then the answer is listed. 


In detail, the program operates as follows. 

Steps 1 § 2 define a vector W of the wavelengths and set up a matrix 
of temperatures . 

Steps 3-7 ask for the input data, multiply it- by X/hc, convert it 
to integer format, and calculate the overall normaliza- 
tion factor, GX. The two GM tube efficiencies are called 
GMl and GM2, although they have been modified. 

Steps 8 S 9 initialize the index I and the vectors for temperature 
(TANS) , blackbody ratios (BANS) , and free-free ratios 
(FANS) . 

Steps 10 - 12 define a matrix X = exp hc/XkT and the free-free spectra 
matrix Y = (W"^) • X. 

Step 13 calculates the free-free ratios corresponding to the tem- 

peratures in T. Step 13 is equivalent to 


FF 


GX 
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Step 14 

Step 15 
Step 16 


calculates 


the blackbody spectra matrix Y, where 


Y = 


X-1 


calculates the blackbody ratios for the temperatures in T. 

increments I, tests it, and cycles the program back to 
step 9, if appropriate. 


Step 17 reshapes the answers for TANS, BANS, and FANS as a matrix 

with the temperature, corresponding blackbody ratio, and 
corresponding free-free ratio appearing in a row; it then 




prints 

them out. 

Step 

18 

prints 

a header for the matrix. 

Step 

19 

clears 

variables . 
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STEPS 1-9 


STEPS 10-12 


STEPS 13-15 


STEP 16 


STEPS 17-19 


Figure F-1 . Flowchart for RATIOS Program 
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V RATIOS 

[1] .25x180 

[2] TEMP^ 3 6 p(0.4 0.6. 0.8 1 1.5 
2 3 4 6 8 10 15 20 30 40 60 80 
100 ) 

[3] CENTER GMli<AREAit( CONVERSION FACTOR):^ 

L4] (?Afl*^L0.5 + (1000 00000)x(JAfli<?jir-«-r/<?ii/l-«-)7x(lf 

1.9 8 62F"8)xG//l^-D 

[53 ‘ ENTER GM2 ^AREA x ( CONVERSION FA CTOR ) ; * ' 

[6] ffM2t-L0.5 + ( 10000000 0)xGM2xJ'fr/GM2't-{/x(li 

1.9 86 2ff“8)x(?A/2-*-a 
[73 GX'^-GXiX 

[83 FANS^BANS^xO 
[93 1^1 

[103 STARTRliI^{li^Z,B^N)o ,^TEMPlIil 
[113 ;?-«-*( ((y<i7o)xj)+(j&i7o)xi'7o) ■ 

[123 2)«.x6pl)^Z 

[13 3 FANS^FANS , X^GX >^ ( GMl + . xjf ) ^ ( (?iV2 + . xj ) 

[143 5)0 .x6pl)4-(Z-l) 

[15] BANS-^BANS , Y^GX^ ( GMl + . xj ) ^ ( GM2 + . xj ) 

[163 -^i^^M-i-D/STARTRl 

[173 U-^RESULT^^{Z 18 p ( ( ,2’SMPxi OOOOOO ) ,5>li75,Fyii75 ) ) 
[183 'TEMP,,BLACKBODY> AND FREE -FREE RATIOS ^ GMl/ GMl . 
[193 TEMP-<rBA NS-^FA NS<r 1 0 

V 
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